#ifndef __SUBLIST_H__
#define __SUBLIST_H__

#include "myIO.h"

inline float Simple_Calc(const ForceInfo& fInfo, const float Meff, const float dist, const float Tm) {
    const float Tm2 = Tm * Tm;
    const float force = -IFr::Potential(fInfo, dist, 'F');
    const float dist1 = dist + force / Meff * Tm2 / 2;
    const float f1 = -IFr::Potential(fInfo, dist1, 'F');

    const float dist2 = dist + (force + f1) / Meff * Tm2 / 2;
    const float f3 = -IFr::Potential(fInfo, dist2, 'F');

    return (force + 2 * f1 + f3) / 4 / Meff;//(force+f1)/2/Meff;
}

#endif // __SUBLIST_H__
